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1 Introduction 

The search for new physics phenomena at the TeV scale begins at the Large Hadron Cohider 
(LHC) this year. Among the proposed models of new physics, extentions of the Standard 
Model (SM) involving supersymmetry (SUSY) are one of the best motivated. Over the years 
much attention has been focused on the Minimal Supersymmetric Standard Model (MSSM) [1], 
characterized by the minimal content of supersymmetric particles and i?-parity conservation. In 
the MSSM the SUSY partners of the SM particles are always produced in pairs. 

The dominant production processes of sparticles at the LHC are those involving pairs of 
coloured particles, i.e. squarks and gluinos, in the final state [2]. The discovery of squarks 
and gluinos should be possible for masses of up to 2 TeV [3]. Depending on the outcome of 
the experimental searches, predictions for the total rates for these production processes will 
either help to determine the masses of the sparticles |lj or to draw exclusion limits for the mass 
parameters. It is therefore important to know the LHC cross sections for production of squarks 
and gluinos with high theoretical accuracy. In particular, control over effects of higher-order 
radiative corrections leads to a stabilization of the theoretical predictions w.r.t. variation of the 
renormalisation and factorisation scales and, consequently, to a smaller theoretical error. 

The leading-order (LO) total cross sections and the corresponding next-to-leading order 
(NLO) SUSY-QCD corrections are known for all hadroproduction processes of pairs of squarks 
and gluinos [SlEllT]. They have been found to be positive and large, especially for the gluino-pair 
production. As pointed out in [7J, an important part of the contributions to the hadronic cross 
sections comes from the energy region near the partonic production threshold. The threshold 
region is reached when the square of the partonic center-of-mass (c.o.m.) energy, s, approaches 
Am?, where m is the average particle mass in the produced pair. The velocity of the produced 



heavy particles in the partonic c.o.m. system /? = ^^1 — Am? /s is then small, /? ^ 1. In this 
region, two types of corrections dominate: Coulomb corrections, due to exchange of gluons 
between slowly moving massive particles, and soft gluon corrections, due to emission of low 
energy gluons off the coloured initial and final states. The large size of the soft gluon emission 
contributions can be traced down, for the perturbative n-th order correction, to the logarithmic 
terms of the form a" log (/3^) where k = 2n, ... ,0. The effects of the soft gluon emission can 
be taken into account to all orders in perturbation theory by performing resummation of the 
threshold logarithms. Resummed predictions are particularly important for processes with large 
masses in the final states, since then the bulk of the production comes from the threshold region. 
This is exactly the case for the production of the sparticles which are expected to be heavier 
than the SM particles. Additionally, if the partonic subprocesses involve gluons in the initial 
state, the soft gluon effects, and thus the impact of resummation, are expected to be significant 
due to the high colour charge of the gluons. 



In this letter, we report on the calculation of threshold-resummed cross sections for hadropro- 
duction of gluino-gluino (gg) and squark-antisquark (qq) pairs. Among the pair-production pro- 
cesses of coloured sparticles at the LHC, the gg production receives the largest NLO SUSY-QCD 
correction [7], and the NLO if -factor, i^NLOi can reach 2 for gluino mass rrig = 1 TeV. The 
corrections to the qq total cross section can be also sizeable i^NLO — 1-3 for the squark mass, 
rriq = 1 TeV, and are the second largest in a certain range of mass parameters. The hadronic 
production of gg and qq pairs are scattering processes with a non-trivial colour flow structure. 
At the level of next-to-leading logarithms (NLL), resummation requires including contributions 
from soft gluons emitted at wide angles. Such emission is sensitive to the colour flow of the 
underlying hard scattering and the evolution of the colour exchange is governed by the one-loop 
soft anomalous dimension matrix [U [9l \TU[ \TT\ [12]. These matrices have been calculated for 
heavy-quark and dijet production O [131 [3 [H] and the general results for any 2 ^ n QCD pro- 
cess with massless particles in the final state have been derived at one- [12], and two-loops |14j . 
Recently, the two-loop anomalous dimension for the pair production of heavy quarks has been 
also determined in the threshold limit [15] . Here we present the explicit form of the soft anoma- 
lous dimension matrices for partonic subprocesses contributing to gg hadroproduction and apply 
the resummation formalism at the NLL level to evaluate the correction due to soft gluon emission 
for the production of gg and qq pairs at the LHC. 

2 Soft gluon resummation 

In this letter we consider the hadronic production of gluino-pair and squark-antisquark final 
states, /ii/i2 — > gg and /11/12 — > qq, where we sum over the left- and right- chiralities of squarks 
and over the spins of gluinos. At LO, these two processes receive contributions from the following 
partonic channels: 

qq -^ 99, qq -^ 99, 99 -^ 99 

and 

qq -^ qq, qq -^ qq, gg -^ qq, 

respectively. The resummation of the soft gluon corrections is carried out in Mellin-A^ space in 
the variable p = Am? / S with S being the square of the hadronic c.o.m energy and m the mass of 
the final-state particle. In the Mellin space, the moments of the partonic cross section ij — > kl 
are given by 

^ij^kl,N{rr?,^J'F,^J'R) = / dp P^'^ ^ij^klip,m^,fJ'F,fJ'R) , (1) 

Jo 

with p = Am? /s. 

The evolution of colour exchange in non-collinear soft gluon emission, which has to be taken 
into account at the NLL accuracy, is governed by the soft anomalous dimension matrix Y^^ [HI 



[9lll0 t[mil2j . In an orthogonal basis in the colour space for which the matrix Y"^^ is diagonal, 
the NLL resummed cross section in the A^-space has the form [9l [13] 



a 



(res) _ ■^'.(0) Pi »i AJ aM) 



"ij-,kl,I,N '^ij~*kl,I ^N+1 ^N+1 ^ij^kl,I,N+l ' V^i 

where we suppress explicit dependence on the scales. The index / in Eq. ([2]) distinguishes be- 
tween contributions from different colour channels. The colour-channel-dependent contributions 
to the LO partonic cross sections in Mellin-moment space are denoted by (^lj_^j^i j ^ and will be 



presented elsewhere |16] . The radiative factors A^ describe the effect of the soft gluon radiation 
collinear to the initial state partons and are universal. Large-angle soft gluon emission is ac- 
counted for by the factors A^'"^^; j ^ which depend on the partonic process under consideration 
and the colour configuration of the participating particles. The functions Cij_kij are in general 
obtained by comparison of the resummed formulas with the fixed-order expressions. They consist 
of Coulomb corrections and A^-independent terms (often referred to in the literature under the 
common symbol Cij^kl)^ containing hard contributions from virtual corrections. For the most 
complete treatment of the threshold effects, the Coulomb corrections should be resummed, see 
e.g. [17^ I13j. The concurrent treatment of the Coulomb and soft gluon corrections is, however, 
beyond the scope of the work reported on here. 

The expressions for the radiative factors in the MS factorisation scheme read (see e.g. 






InAV = / dz"— I ^Ai{as{q'^)) . 



In Agl,^,^^ = j^ dz -^—— A,.fc.,/(as(4m2(l - zf)) . (3) 

The coefficients T = Ai, Dij^f^ij are power series in the coupling constant ag , J-' = {^)J-'^^^ + 
(^!^^2jr{2) _^ _ __ rj.Yie universal LL and NLL coefficients a\^\ ^f ^ are well known [iHl [19] and 
given by 

with Cg = Ca = 3, and Cq = Cp = 4/3. The customary NLL expansions of the radiative factors 
can be found in [13J. Since after the expansion of the exponentials the terms constant in N 
contained in Cij^ki generate contributions of the next-to-next-to-leading logarithmic (NNLL) 
order, and we do not include Coulomb corrections here, we keep Cij_kij = 1 for the purpose of 
the calculations. 

The NLL coefficients D . -^ , are the same as in the heavy-quark production and were 
calculated in [9l [13] 

D^^^ - = D^^^ - = 

D^^} -, = D^^^ -, = -Ca, (5) 

gg^m,2 gg^qq,^ ^' ^ ^ 



with the index 1 = 1,2 corresponding to the singlet and octet {1,8} exchange in the s-channel. 
For the gg production we first study the soft anomalous dimension matrices V^^^ds and then de- 
rive the coefficients D\-_^~~ r. Apart from obvious colour-dependent modifications, the matrices 
Y'^3^99 are defined in the same way as the soft anomalous dimension matrices for the heavy- 
quark production process in [^. For the purpose of this calculation we consider the partonic 
processes 

q{pi,ai) q{p2, 02) ^ g{P3, 03) 9{Pi, «4) 
and 

g{pi,ai) 9{P2,a2) -^ g{p3, 03) g{p4, 04) , 

where pi are the particle four-momenta and Oj and at are the colour indices in the fundamental 
and adjoint representation of SU(3), respectively. 

Three independent tensors are forming the basis in the space of colour exchanges for the 
QQ ~^ 99 process. We choose an orthogonal s-channel basis {cj}, I = 1,2, 3: 

where T matrices are the SU(3) generators, and the tensors correspond to 1, 8s and 8a 
representations, respectively. For the gg channel there are eight independent colour tensors. 
Following [TTj we choose an orthogonal basis {c^} (I = 1, 2, . . . , 8) consisting of five tensors cf, 
cf, C3, C4, Cg corresponding to the 1, 8s, 8a, 10 © 10, 27 representations in the s-channel, and 
three additional tensors Cg, Cy, Cg. The base tensors are 

8 

Q 

9 _ ^^aia2bjba3a4 

2 5 

„9 _ faia2b rba^a^ 

^9 _ _ r^aia3^a2a4 _ ^aia4^a2a3\ _ ^ taia2b ^basa^ 

1 1 3 

fS — Z. C^aia3^a2a-4 _|_ ^aiai ^a2a3\ _ _^aia2^a-ia4, _ ^^01026^60304 

2 8 5 
g9 _ 1 / ^1026^60304 _|_ ^01026^60304 

g5 _ 1 ^"^1026^60304 _ ^01026^60304 

g9 _ 1 ^^01036 ^60204 _|_ ^oi036^6a204\ /j\ 



We introduce the notation 

T = In 

C7 = In 



m^ — t\ 1 — in 



^ \Jw?s ) 2 

rr? — u\ 1 — iiT 

\rrrfl't 



s 



Lp + l 



S ^ -^^ , (8) 

4 



where the Mandelstam variables are given by 



(Pl+P2)^ i={pi 



P3) 



u= {pi- piY 



and Lp = ^(1 - 2m^/s) (in j^ + iii). We also define A = T + U , Vt = T-U . 

Using the results for the one-loop eikonal integrals from [9] we calculate the matrices V^^^dg JJg 
For the qq channel we obtain 



Yii-'99 



vr 



/ 



V 



6S 



-2f] 





35 + |A 



-n 

35 + |A 



\ 



-zvrl 



(9) 



The matrix for the gg channel has the block form 



rv 

Y99-ag — zl 

■K 



where the five-dimensional matrix reads 

/ 




3i7rl 



(10) 
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35 + f A 





!^ 










3A 


19 
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1^ 


^9 


4A - 25 y 



(11) 



and the three-dimensional matrix Fs is diagonal, 

Fs = diag(3(5 + C7), 3(5 + f), 3(f + C7) ). 



(12) 



As in [9j, the matrices F contain terms from the one-loop eikonal integrals for 2^2 scattering 
with massive particles in the final state and are shifted by half of the soft anomalous dimension 
for the Drell-Yan cross section Q|. The matrices Q , (|10p complement the results of [TT] for the 
case of pair-production of massive colour-octet particles with equal masses. 

At the threshold, s — > Arm?', the soft anomalous dimension matrices approach the diagonal 
form, 

/y 

p,,-^ss^^diag(7f,7f, •••,7f), 
vr 

^x '■ ' q q q\ 



n'?-'?'?^_idiag(7?,72^7l) 



(13) 



^ Gauge dependence in the one-loop integrals cancels against gauge dependence of the incoming-jet factors 
(which corresponds to gauge dependence of the soft function in the resummed Drell-Yan cross section [111 123] ). 
leaving the results presented here gauge-invariant. 



The off-diagonal terms, proportional to 0,, vanish like f3 for /? ^ and may be neglected. Using 
^99^9,1 = 2I^«(^/) ^"^d ^S^99,i = 2Re(7f) [inilll] we obtain 

{^^^9, /} = {0, -3, -3, -6, -8; -3, -3, -6} 

{<-U/} = iO'-3'-3}- (14) 

Note that the values of the D'-^-'-coefficients are the negative values of the quadratic Casimir 
operators for the SU(3) representations for the outgoing state. This agrees with the physical 
picture of the soft gluon radiation from the total colour charge of the heavy-particle pair produced 
at threshold [13]. 



3 Numerical results 

In the phenomenological analysis we consider a wide range of gluino and squark masses. Left- and 
right-handed squarks of all flavours are assumed to be mass degenerate. For the gg production 
we varyn the gluino mass, rrig, between 200 GeV and 2 TeV. Similarly, for the qq production we 
take 200 GeV < rriq < 2 TeV. We present the results for a fixed ratio of gluino and squark masses, 
r = — 2-, and choose the following values r = 0.5, 0.8 1.2 1.6, 2.0. The qq cross section accounts 

Trig 

for production of all qq flavour combinations apart from the ones with scalar top particles. 

The main phenomenological results of this letter are the resummation-improved predictions 
for the pp -^ qq and pp -^ gg total cross sections at \/S = 14 TeV. The resummation-improved 
cross sections are obtained through matching the NLL resummed expressions with the full NLO 
cross sections 

rCMP+ioo ^jy 
^CMP—ioo 



("^"*"^) (n 77,2 r 2.^ _ ^^ r^MP+^oo ^ f (^+1)^,^2 ^ .(^+1)^,^2 ^ 

■ ' q,q,g ''Cmp-ioo ^TTt 



«J=9,'?,9 



%- 



(NLO) 



.(NLO) t^ _2 J, .2 



+ <riH(P,-M/«^}), (15) 

where {fi^} = {np, fij^}, cr^jZlij^ is given in Eq. @ and (T-j^^^^^ \(i<ilo) represents its perturbative 
expansion truncated at NLO. The moments of the parton distribution functions fi/h{x, l^-'p) are 
defined in the standard way 



^ u 



/o 

The inverse Mellin transform (|15p is evaluated numerically using a contour in the complex-A^ 
space according to the "Minimal Prescription" method developed in Ref. [20j . 



^For the highest masses considered here, the experimental exploration will require the integrated luminosity 
of 0(100 fb-^). 



0.2 



0.15 



0.1 



0.05 



r = 0.5 




r = 0.8 




r = 1.2 




J 


r = 1.6 


/ 


// 


r = 2.0 


j/^\ 




/^^ •'' 




// .•' ' 




//. >> 




// .-.'y 




y/..'y 




X'.:'y 




•O^v'.^ 




yAi-y 


. 


<>uy 


>/-' -•'^ 1 




V 


y^'^ 


. 


y^ 






a) 



0.04 



0.03 



0.02 



0.01 



r = 0.5 


■ 




/ 






r = 0.8 


— 


-- 


/ .■• 


r = 1.2 
r = 1.6 
r = 2.0 






/ /yy 






— 


~ / 






r 








b) 



0.2 0.5 1 1.5 

mg [ TeV ] 



0.2 0.5 1 1.5 

mq [ TeV ] 



Figure 1: The relative NLL correction i^NLL — 1 for the gg (a) and the qq (h) total production 
cross section at the LHC as a function of gluino and squark mass, respectively; r = rrig/mq. 

The NLO cross sections are evaluated using Prospino [21j, the numerical package based 
on calculations employing the MS renormalisation and factorisation schemes. We use the 
CTEQ6M [22] parameterization of parton distribution functions (pdfs) for all numerical pre- 
dictions. Similarly to other available pdfs, the assumption of five massless quark flavours active 
at large scales is made in the CTEQ6M parameterization. Consequently, in the NLO and 
NLL calculations we use the two-loop MS QCD running coupling constant Os with rif = 5 and 
Aqqj-, = 0.226 GeV. The effects due to virtual top quarks and virtual sparticles in the run- 
ning of Os and in the evolution of pdfs are thus not included in our predictions. However, the 
value of the top mass, rrit = 175 GeV, enters the matched NLL cross sections through the NLO 
corrections. 

All numerical calculations were performed using two independent computer codes. As a first 
check we reproduced numerically, and analytically, the LO results OET]. We also verified that 
the NLL resummed expression in A^-space generates the same logarithmic terms as these present 
in the Mellin transforms of the NLO correction in the threshold limit given in [7] . 

The relative corrections from soft gluon resummation to the gg and qq NLO production 
cross sections at the LHC, 
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Figure 2: S'ca/e dependence of the total gg (a) and qq (b) production cross section at the LHC 
(see the text for explanation). 



are presented in Fig. la and Fig. lb, respectively. In the plots we set the scales fip = ^_R = MOi 



where ^o = "ig (^o 



rrifj 



for the gg production (the qq production). -STnll grows with the 



final-state mass and depends on the mass ratio r in a moderate way. The correction, i^NLL — Ij 
reaches 16% (8%) for the gg production with r = 1.2 and nig = 2 TeV (1 TeV), and 4% (2%) 
for the qq production with r = 2 and rUg = 2 TeV (1 TeV). The stronger effect found in the 
gg production follows from the dominance of the gg — > gg channel, and hence more intense soft 
gluon radiation. 

We also investigate the dependence of the matched NLL cross section on the values of 
factorisation and renormalisation scales, in comparison to the NLO cross section. To illustrate 
our results we choose n = ^p = l^R and r = 1.2. In Fig. 2a and Fig. 2b we plot the ratios 
^NLO(^ = i^iQ)/a^^^{ix = /io) and a('^^*=^^) (^u = i^lQ)/(J^'^''^''^\^l = /xq), obtained by varying ^ 
between ^ = 1/2 and ^ = 2, for gg (^uq = mg) and qq (^o = mtq) production, respectively. Due 
to resummation, the scale sensitivity of the gg production cross section reduces significantly, by 
a factor of ~ 3 (~ 2) at rUg = 2 TeV [rUg = 1 TeV). At rUg > 1 TeV the theoretical error of 
the matched NLL gg cross section, defined by changing the scale ^ = ^p = IJ'R around hq = rUg 
by a factor of 2, is around 5%. In the case of the ijij production, the reduction of the scale 
dependence due to including soft gluon corrections in the theoretical predictions is moderate. 
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